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MATHEMATICAL DETERMINATION 
OF GEOMETRICAL IMAGE ABERRATIONS IN 
SINGLE- AND DOUBLE-MIRROR SYSTEMS 

By Martin G. Hurwitz and Hubert F. A. Tschunko 
Electronics Research Center 


SUMMARY 

The mathematical analyses to determine the image aberrations 
of single- and double-mirror systems are presented. The single 
mirrors, spherical and paraboloidal, are considered. Then, 
double-mirror systems with paraboloidal primary and hyperboloidal 
secondary are considered, along with higher aspherics of the 
secondary and tolerances in the relative positions of the mirrors. 
This investigation is a first step in the evaluation of the 
imaging performance of single- and double-mirror systems. 

The general analysis available in the literature is not 
applicable without development for the special requirements of 
the present case. This special raytracing is the prerequisite 
needed later for the Fourier transformations which serve to derive 
the image performance data necessary to the determination of wave 
optical aberrations. 


INTRODUCTION 

Image aberrations of mirrors and mirror systems have been 
investigated by many authors. In general, the more recent in- 
vestigations used the third and higher order theory of aberra- 
tions. In this treatment these theories are not used, and the 
different families of geometrical image aberrations are not con- 
sidered . 

The image aberrations of single- and double-mirror systems 
are analyzed. Three-dimensional analytical and numerical ray- 
tracing methods are used in the determination of the sizes and 
shapes of the images. Apertures are represented by points uni- 
formly distributed on concentric aperture circles which are pro- 
jected as a cluster of closed-image curves. The largest dimen- 
sion of each cluster is taken to be the aberration value e. 

First, single imaging mirrors with both spherical and 
paraboloidal imaging surfaces are investigated. Special features 
of the image curves, including the catacaustic for the sphere, 
are considered. Then, double-mirror systems with paraboloidal 
primary and hyperboloidal secondary mirrors are treated. Finally, 
position tolerances in the distances between both mirrors are 
considered. 



MATHEMATICAL ANALYSIS 


Single Mirrors 

Initially, single -mirror systems- were considered using 
mathematical models of the spherical and paraboloidal mirrors. 
These models were developed to permit comparative investigations 
of aberrations of each mirror as a function of its focal number 
and angle of field. 

Consider a three-dimensional Cartesian reference frame with 
the optical axis along the z axis. A single-mirror telescope is 
modeled by a portion of a sphere or paraboloid with its vertex 
at the origin, its axis coincident with the optical axis, and its 
center of curvature along the positive z axis (Figures 1 and 2) . 
The equations of the surfaces are 



Figure 1 Figure 2 

A point target at infinite distance emits light rays imping- 
ing on the aperture; these rays are reflected from the mirror 


2 



y 



surface and imaged at a focal 
plane placed orthogonally to the 
optical axis at the focus of the 
paraboloid (0, 0, 0.5) where the 
focal length F L is 0.5. The 
apertures, in either case, are 
represented by points uniformly 
distributed on concentric aper- 
ture circles representing differ- 
ent focal numbers (Figure 3) . 

In all the ray tracing cal- 
culations, directions are defined 
in terms of 


Figure 3 
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and cos y. 


These are the direction cosines 
related to the x, y, and z axes, respectively, and the subscripts 
j = i, ; n, or e are related to incident, normal, and exit rays, 
respectively. The exit ray directions may be determined when the 
incident rays and surface normals are known. Thus: 
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Consider a ray approaching the mirror with angles y 
= tt/2 , = (tt/ 2) - y^. 

The intersection of the mirror surface with its normal is 
determined from the focal number, i.e., focal length over aper- 
ture diameter equals F N , as follows: 


y 0 


0.5 


2F 


N 


(4) 


where y 0 is the radius of the aperture. The radius of curvature 
is assumed to be 1.0 and the focal length equal to 0.5. For the 
sphere : 


*0 


sin 


1 n 


(5) 
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and for the paraboloid: 


y 0 = _tan Y n - < 6 ) 

The direction cosines of a , B are determined from 

n n 

cos a n = sin Y n cos ip , cos 8 n = sin Y n sin ip (7) 

where i p ranges as position angles over the aperture circle points. 
The direction cosines of the exit ray may now be determined from 
Eg. (3). 


The exit ray is defined in terms of its intersection with 
the mirror (x^, y^, z^) , and its direction cosines as 


where 


while 


x - x 1 y - y-L z - z 1 

cos a cos 6 cos y 

e e e 


(8) 


x i = y o cos ^ ' Y 1 = y o sin ^ 


(9) 


z 1 = 1 - cos y n 


(10) 


for the sphere and 


If 2 ^ 2\ 

: 1 2 \ X 1 + y lj 


(ID 


for the paraboloid. 

This ray may be intersected with the focal plane z = 0.5 
to define an image point. The entire image is obtained by 
allowing tj; to range over each aperture circle of constant radius. 
Different focal numbers provide different radii and different 
concentric aperture circles represent different focal numbers. 

All image curves together represent the image cluster, the 
whole geometrical image of the point object. 
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Special Image Points 


The images from a paraboloidal single mirror display some 
interesting properties, such as a triple point and highest point, 
already pointed out by Crockett (ref. 1) and Plummer (ref. 2). 

The triple points indicate a region in focal numbers and 
fields of view where the maximum sizes of the geometrical aber- 
rations change from linear to quadratic dependence with field of 
view. 

Each highest point indicates one limiting position within 
an image cluster. These points are necessary in determining the 
aberration sizes of the clusters, especially for larger focal 
numbers. Without these points no theoretical evaluation of the 
aberrations for higher focal numbers is possible. 

These derivations are included because no exact determina- 
tions of these points may be found in the literature. A deriva- 
tion of these points deals with relations between the angles of 
the incident ray and normal rays and is presented below. 

Triple Point — The triple point is that point at which the 
image curve intersects itself for \p = 0° and ±90°. 

Let \p = 0°. Then, from Eq. (7) : 

cos 6 =0, cos a = sin y (12) 

n n n 

and 

cos y^ = 0, cos = -sin y^. (13) 

From Eq. (3) : 
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The intersection points (x^, y.^, z^) of the incident ray 
with the mirror surface are found to be 

1 2 

x x = -tan y n , y 1 = 0 , z ± = ^tan r n (15) 
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and the exit ray is then 


x - x x y - Y-l _ z - z ± 

cos a cos 8„ - cos y^ 

e e e 

Thus, at z = ■— : 

cos a / \ 

x = ssst ;! 2 - z i) + x i 

- tan 2 v n(l - ! tan2 Y n) - tan T n 
= 0 


( 16 ) 
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Now, let ip = -9 0 . Then: 


cos 3 n = -sin Y n , cos a n = 0 


cos a. = 0, cos 8* = -sin y. 


From Eq. (3) : 


cos a =0, cos 8 = sin (y. - 2y ) 


n 


cos y s = cos (y^ - 2y n ) . 


(17) 


(18) 
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The intersection of the incident ray with the surface now becomes 


x 1 = 0 


y, = -tan y 
J 1 'n 


At 
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( 20 ) 


The two cases, ip = 0 and ip = -90°, both yield x = 0. The 
two values of ip result in the same image point, so the y values 
must be equated: 
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( 21 ) 


This is simplified as follows: 


(tan YjJ (tan 2y n ) (tan - tan y n ) = 0 (22) 

Thus, the possible solutions: 

tan Y i = 0, tan 2y n = 0, tan y^ = tan y n (23) 
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yield the values 


y ^ = ±mr 


Y 


n 



n = 0, 1, 2 , . . . 


n = 0, 


1 , 2 , . . . 


(24) 


and 


tan y ■ 
_x 

tan y n 


(25) 


Equation (25) provides the condition for the occurrence of the 
triple point. 

Highest Point — The highest point of the image curve of a 
paraboloid mirror is found as follows. 

Let the highest y coordinate of the image of a ray be 
denoted by y* and the corresponding point of intersection of the 
ray with the mirror by y^. The highest point is then determined 
by minimizing y* - y. 

From Figure 4: 


tan Y e = tan (2 Yn - y ± ) 


Ay 2 + 2y - A 
-y 2 + 2Ay + 1 


(26) 


where 


y = tan y 
2 ' n 


and 


(27) 


A = tan y ^ • 
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The value: 


Figure 4 


y* - y = tan Y e (z p - z x ) , 


(28) 


with z F equal to the position of the focal plane, 
minimized . 

Taking 



is to be 


and tan y as above , 
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Y* - y 
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(29) 
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For a minimum: 


-A(y* - y) = 0. (30) 

This results in the equation: 

y 5 - 3 Ay 4 - 2y 3 - 2Ay 2 - 3y + A = 0 (31) 

with 


y / A ± \/a 2 + 1 . (32) 


Then A can be found as 


tan y ± 


A = Y 5 ~ 2y 3 - 3y 
3y 4 + 2y 2 - 1 


(33) 


When y = tan Y n is small, A « 3y, 
or 

tan y . « 3tan y 
' i ' n 


or 

tan y i 

lil * I 

tan y^ 3 


(34) 


(35) 


for small y , and deviates from 1/3 as y increases. 

'n' ' 'n 

With this result, the highest point y* may be determined 
from the above. This equation (35) is the condition for the 
occurrence of the highest point. 

Minimal Image Di a meter of the Reflected Ray Bundle in the 
Catacaustic of a Spherical Mirror — In the case of the spherical 
single mirror, the incident rays form a blurred image created by 
the intersections of the tangents to the catacaustic curve with 
the optical axis, as shown in Figure 5. 
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The catacaustic may be described parametrically as: 


z = £-(3cos a - cos 3a) 

4 

(36) 

y = j(3sin a - sin 3a) 

These equations may be combined to yield 

[(Mf + (^f- 4 ] = 1°8(^) 2 (37) 

where r is the radius of the sphere, and will be taken to be 
unity. 

Further reduction of the above yields 


y 


2 



(38) 


The reflected rays from a spherical mirror create a blurred 
image even when the incident rays are parallel to the optical 
axis. This occurs because the sphere does not image to a single 
focus as does the paraboloid, but each concentric zone of the 
spherical mirror has a different focal point for angles of 
incidence yi = 0. In this way, the catacaustic is formed as an 
envelope of the reflected rays. The rays from the largest 
aperture will be imaged in the focal plane, with the largest 
aberrations. The reflected ray is given by 


y = -tan Y e (z ~ Z]_) + Y]_ 


(39) 


For the sphere: 


Y = 2 y 
'e n 


y-L = sin Y n 


(40) 


z 1 = -cos y n 
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Substituting : 


y = -tan 2y n (z + cos y n ) + sin Y n • 


Solving for z: 


y - sin Y 
2 n 

-tan 2y 
' n 


cos y 


n 


(41) 


(42) 


Substitution of this expression for z into the equation of the 
catacaustic then yields: 


E m 1 - 0 

i=0 

where 


R 0 - (« 2 - 

*1 ■ ¥( b 2 - kf 


, 3 / 

f _ 4 

3B 2 

1 \ 


/ 4 B 2 

5 \ 

^2 2 

, 5B - 

2 

■ + Is) 

+ 3 

( B - T - 

64) 

A = -L| 

3 t 3 

(20B 2 

- 3 ) 

+ 12B ( 
T ' 

B 2 

-\) 


3 

/ 2 

i\ 


2 

1 \ , - /_2 

1 \ 

A . = —r- 

4 rp4 

\5B - 

?) 

12 ( 6E 
T 1 

r - 

2) + 3 ( B 

4 ) 
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(43) 
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and 


B 


-sin Y + cos y tan 2y 
' n n 'n 


tan 2 y 


n 


(45) 


T = tan 2 y 

n 


(46) 


This equation in the variable y is solved by an iterative 
procedure, while noting that y and yn are opposite in sign. The 
diameter of the minimum bundle is taken as the aberration. 


e = 2 | y 


(47) 


To this minimal bundle diameter belongs a z-coordinate 
defining the optimal focal position for Yi = 0- This focal 
position yields the focal length used to determine the image 
aberrations of the spherical mirror. For each focal number F N 
a nominal value of the largest aperture diameter is used, thus 
yielding the optimal focal position length, or optimal focal 
length. For all smaller aperture circles of the same aperture, 
the same optimal focal length had to be used to assure the same 
focal plane in the same cluster of image curves. 

The optimal focal plane positions OFP are presented as a 
function of focal number in the following table. 

The OFP ' s are often obtained in a cumbersome and rather 
imprecise way. For this reason, they were computed with a higher 
number of figures to show the approach to the value of 0.5. 

These values are necessary to compute the image aberrations 
of spherical mirrors over a larger range of focal numbers. 


DOUBLE-MIRROR SYSTEM 

First the ideal system is treated. This is a system of 
geometrically perfect paraboloidal primary and hyperboloidal 
secondary mirrors, without any higher surface aspherics, nor any 
provisions for tolerances in the mutual position between the two 
mirrors . 
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VALUES OF OPTIMAL FOCAL PLANE AS A FUNCTION 
OF FOCAL NUMBER FOR SPHERICAL MIRRORS 


F 

N 

OFP 

F 

N 

OFP 

.5 

0.43193 40028 65542 

25.0 

0.49998 12482 41999 

.75 

0.47669 87802 44569 

32.0 

0.49998 85552 53325 

1.0 

0.48754 48957 61076 

50.0 

0.49999 53123 90134 

1.5 

0.49465 18525 90037 

64.0 

0.49999 71389 36123 

2.0 

0.49702 66638 91014 

80.0 

0.49999 81689 28549 

2.5 

0.49810 72309 29950 

100.0 

0.49999 88281 18134 

3.0 

0.49868 93758 44114 

150.0 

0.49999 94791 65315 

3.333 

0.49879 90682 39857 

200.0 

0.49999 97070 30821 

4.0 

0.49926 48846 15958 

250.0 

0.49999 98124 99825 

5.0 

0.49953 01484 09819 

256.0 

0.49999 98211 86497 

6.0 

0.49967 39483 58068 

320.0 

0.49999 98855 59363 

8.0 

0.49981 67267 17227 

400.0 

0.49999 99267 57878 

10.0 

0.49988 27437 89342 

500.0 

0.49999 99531 34400 

12.0 

0.49991 85866 62538 

512.0 

0.49999 99552 96535 

16.0 

0.49995 42131 52687 

1000.0 

0.49999 99884 33952 

20.0 

0.49997 06988 32746 

1024.0 

0.49999 99889 19132 


Ideal System 

The double-mirror systems are Cassegrain systems. The 
primary is a paraboloid of revolution and the secondary is one 
sheet of a hyperboloid of two sheets. The paraboloid contains a 
central core in an area obstructed by the projection of the 
secondary mirror in the direction of Yi, and of sufficient 
diameter to permit imaging of rays in a focal plane described 
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by z = -B, B > 0. Except for the core, the mirror will be the 
same as the single-mirror paraboloid described above. 

In the ideal system, one focus of the hyperboloid coincides 
with that of the paraboloid, while the other focus is located at 
the point (0, 0, -B) , the intersection of the focal plane with 
the optical axis (Figure 6) . 

The hyperboloid is described by the equation 


z = h + j[ (K + 1) (z - h) 2 + p 2 ] (48) 

where 

p 2 = (x - x Q ) 2 + (y - y Q ) 2 , (49) 

x Q and Yq are displacements from the optical axis. 


C 


2v c K 
(K + 1 ) L 


(50) 


is the vertex curvature of the hyperboloid, 

T , /M + 1\ 2 _2 _ . , „ 

■K — 1/ ® ® > 1 , (51) 

where e is the eccentricity of the hyperboloid, and M is the 
imaging magnification of the secondary mirror, 

L = 0.5 + B , (52) 


which is the distance from the secondary effective focal plane 
to the primary focus, and, 


0.5 - MB 
M + 1 ' 


t 

(53) 


is a displacement along the optical axis. 
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PRIMARY MIRROR 
/ F n = 4 

POSITION OF 



I 


The hyperboloid equation may be written as 



where 


a 


1 

C(K + 1) 


( 54 ) 


(55) 

r 


b = a^- (k + 1) 


(56) 

t 


are the semi -major and semi-minor axes, and 


z q = h + a (57) 

is the z coordinate of the center of the system. This provides 

a more common form for an elliptic hyperboloid of two sheets with 

center at (x , y , z ) . 

o' *o' o 

A further simplification yields: 


z = 1 V * 2 - kTT + z o (58 > 

with the plus sign representing the right-hand sheet. 


Double Mirror Geometry 

The double-mirror system will create images in the back 
focal plane. The imaging process will be initially the same as 
that for the single mirrors up to reflection from the primary. 

The exit rays from the primary will intersect the secondary at 
points which must be determined. This determination is accom- 
plished by simultaneous solution of the equation of exit ray with 
that of the hyperboloid using an iterative procedure. The inter- 
section point will be called (x 2 , y 2 , z 2 ) , and the value for z 2 
should approximate the value B + h. 

Once this point is established, the ray is reflected from 
this point, through the core of the paraboloid, and toward the 
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focal plane. The normal to the hyperboloid, H(x,y,z) = 0, at the 
intersection point is established for use in determination of the 
exit ray. This normal will have direction cosines 


3H. 

cos “n = 9^ 

cos S = l^t 
n 8 y 

COS Y = -^-S-t 

cos Y n 3z u j 


(59) 


with t determined as 


t = 


Vdf) 2 ♦ (n 1 ♦ (H) 2 


(60) 


since 


2 2 2 
cos a + cos 6 + cos y = 1 . 

n n 1 n 


(61) 


The equation of the exit ray from the hyperboloid is then 


x - x 2 _ y - y 2 _ z - z 2 

cos a cos 6 cos y 
e e ' e 


(62) 


where cos a e , cos 6 e , and cos y e are determined as in Eqs . (3), 
using the old exit angles as the new incidence angles, and the 
normals just established. The intersection point of this ray 
with the focal plane z = -B is then established, providing an 
image point (x^, y 3 , Z 3 ) . 

The values y 2 may be reduced to a common origin so that 
comparison of values is facilitated. This is accomplished by 
subtracting a normalizing value from each y 3 . Such normalizing 
values are computed as above with 

Y n = °/ Y ± 7* 0, and F N = F N (nominal) 
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and represents the intersection point of a ray through the center 
of the aperture with the secondary effective focal plane. 

When this procedure has been done for all values of ip for 
each aperture, the required images have been achieved as clusters 
of closed image curves. The largest dimention of each cluster 
is taken to be the aberration value e. 


Higher Aspherics 

The hyperboloid will be described so as to include higher 
aspherics. These higher aspherics are introduced in such a way 
as to preserve radial symmetry of the surface about the optical 
axis. The inclusion of higher aspherics requires a modification 
of the equation of the surface, 


z 



2 

P 


a 2 (K + 


1 ) 


n 

+ Z 

i=l 


c iP 


2 (i + 1) 


(63) 


with the C^'s as parameters to be determined. For this investi- 
gation, n is taken to be no more than 4 because the ranges of 
values to be considered render negligible the terms p 1 with 
i > 10. 


Tolerances 

In addition to higher aspherics, the hyperboloid may be 
repositioned. A position change along the optical axis, an axial 
shift, will be denoted by 6, and the hyperboloid equation becomes 


z = Zq + 6 + a-t/1 


4 


n 


a (K + 1) i=1 


£ 


2 (i + 1) 


(64) 


and represents the secondary mirror, an aspheric hyperboloid dis- 
placed axially by an amount 6 and laterally by amounts x q and y Q , 
where x and y are now to be termed lateral shifts from the 
optical°axis , as defined in Eq. (49). 
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The non-displaced case without higher aspherics is given by 
Eq. (64) with 6 = x Q = y Q = 0, and = 0, i = 1,...,4. 

Finally, the hyperboloid may be tilted about its vertex. 
This tilting will be considered relative to the x and y axes, 
respectively, and jointly. Two-dimensional Cartesian coordinate 
systems, with x, y, and x" , y' axes rotated through an angle 8 
relative to each other may be related as 


x = x'cos 8 + y'sin 8 

y = -x'sin 0 + y'cos 8 . (65) 


To achieve rotation of the system, the on-axis hyperboloid is 
shifted along the optical axis a distance, 


A 


h + 


2 

C (K + 1) 


+ 6 


( 66 ) 


so that its origin coincides with that of the coordinate system. 
The hyperboloid is rotated first about y with angle 8 , then 
about x with angle <j> , and finally shifted back to its original 
vertex position. The equation of the hyperboloid with tilting 
is 


[ (z - A)cos <j> + (x - x Q )sin <j>]cos 0 + (y - y Q )sin 8 


= a 




+ Ip-UL 


->] 


+ f (p* ) 


(67) 


where 


n 


f(p') = Z c i<p’> 2(i + 1) 

i=l 


( 68 ) 
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and 


n O 

(p 1 ) = [(^ - z)sin + (x — x q )cos 

+ |-[(z - A)cos <t> + (x - x Q )sin <f>]sin 6 
+ (y - y Q ) cos e } 2 (69) 

Equation (67) will be rewritten as 

H (x,y , z) = 0 (70) 


where 


H(x,y,z) = [(z - A)cos 43 + (x - x Q )sin <j)]cos 0 


+ (y - y ) sin 6 - a 




(P ') 2 


- 1 


- f(p') 


(71) 


Optimum Axial Focus with Axial Shift Tolerance 

In an ideal, double-mirror system, a parallel bundle of 
light approaching the primary mirror with yi = 0 would create a 
point image on the optical axis in the focal plane. (When 

^ 0, aberrations are formed in the focal plane image.) If, 
in addition, tolerances in axial shift are introduced, it becomes 
necessary to find the optimal position for the focal plane, as 
the image with minimal aberrations no longer occurs at Z = -B, 
the secondary focus of the original hyperboloid. 

The optimal focal plane position was evaluated as follows. 
For $ = 90°, the corresponding ray was traced as above, and 
imaged in each of two focal planes separated by a distance Zg. 
This was done as a function of the focal number F N ; the ray 
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corresponding to the minimum value of F N was intersected with 
the rays corresponding to Fjg + nA (Fjg) , with A (F N ) an increment, 
for n = 1,2,... until an intersection was obtained which pro- 
vided the minimum aberration value. 

An iteration was then performed to ensure the realization 
of the actual minimum aberration. The position on the optical 
axis corresponding to these minimum aberrations was taken to be 
the focal plane position. 


COMPUTATIONS 

The given analyses provide the necessary geometrical skele- 
ton for the numerical evaluation of image aberrations for the 
implementation on digital computers . 

The computations were done using FORTRAN IV on the IBM-7049 
II Of NASA-ERC. 
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